Weak measurements are universal 
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It is well known that any projective measurement can be decomposed into a sequence of weak 
measurements, which cause only small changes to the state. Similar constructions for generalized 
measurements, however, have relied on the use of an ancilla system. We show that any generalized 
measurement can be decomposed into a sequence of weak measurements without the use of an 
ancilla, and give an explícit construction for these weak measurements. The measurement procedure 
has the structure of a random walk along a curve in state space, with the measurement ending 
when one of the end points is reached. This shows that any measurement can be generated by 
weak measurements, and hence that weak measurements are universal. This may have important 
applications to the theory of entanglement. 



In the original formulation of measurement in quantum 
mechanics, measurement outcomes are identified with 
a set of orthogonal projection operators, which can be 
thought of as corresponding to the eigenspaces of a Her- 
mitian operator, or observable [H, Q- After a measure- 
ment, the state is projected into one of the subspaces 
with a probability given by the square of the amplitudc 
of the state component in that subspace. 

In recent years a more general notion of measure- 
ment has become common: the generalized or positive- 
operator valued measurement (POVM) This formu- 
lation can include many phenomena not captured by pro- 
jective measurements: detectors with non-unit efhciency, 
measurement outcomes that include additional random- 
ness, measurements that give incomplete information, 
and many others. POVMs have found numerous applica- 
tions, especially in the rapidly-growing field of quantum 
information processing [J|. 

Upon measurement, a system with density matrix p 
undergoes a random transformation 

Pj =M jP Mj/ Pj , ^TaíJaí, =1, (1) 
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with probability pj — Tr(MjpMj), where the index j 
labels the possible outcomes of the measurement. This 
transformation is commonly comprehended as a sponta- 
neous jump, unlike unitary transformations, for example, 
which arc thought of as resulting from continuous unitary 
cvolutions. Any unitary transformation can be implc- 
mented as a sequence of weak (i.e., infinitesimal) unitary 
transformations. One may ask if a similar decomposition 
exists for generalized measurements. This would allow us 
to think of POVMs as resulting from continuous stochas- 
tic evolutions and possibly make use of the powerful tools 
of differential calculus in the study of the transformations 



that a system undergoes upon measurement. 

In this paper we show that any generalized measure- 
ment can be implemented as a sequence of weak measure- 
ments [a, [6( . We call a measurement weak if all outcomes 
result in very small changes to the state. (There are other 
dcfinitions of weak measurements that include the possi- 
bility of large changes to the state with low probability; 
we will not be considering measurements of this type.) 
Therefore, a weak measurement is one whose operators 
can be written as 

M j = q j (ï + e j ), (2) 

where < qj < 1 and è is an operator with small norm 
||e|| -C 1. Weak measurements have been studied both in 
the abstract, and as a means of understanding systems 
with continuous monitoring. In the latter case, we can 
think of the evolution as the limit of a sequence of weak 
measurements, which gives rise to continuous stochastic 
evolutions called quantum trajectòries [g, 0| ■ 

It has been shown that any projective measurement 
can be done as a sequence of weak measurements; and 
by using an additional ancilla system and a joint uni- 
tary transformation, it is possible to do any generalized 
measurement using weak measurements [8]. This pro- 
cedure, however, does not decompose the operation on 
the original system into weak operations, since it uses 
operations acting on a larger Hilbert space — that of the 
system plus the ancilla. If we wish to study the behavior 
of a function — for instance, an entanglement monotone — 
defined on a space of a particular dimension, it compli- 
cates matters to add and remove ancillas. We will show 
that an ancilla is not needed, and give an explicit con- 
struction of the weak measurement operators for any gen- 
eralized measurement that we wish to decompose. 

It is easy to show that a measurement with any number 
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of outcomes can be performed as a sequence of measure- 
ments with two outcomes. Therefore, for simplicity, we 
will restrict our considerations to two-outcome measure- 
ments. To give the idea of the construction, we first show 
how every projective measurement can be implemented 
as a sequence of weak generalized measurements. In this 
case the measurement operators P\ and Pi are orthogo- 
nal projectors whose sum Pi+P 2 = J is the identity. We 
introduce the operators 

P{x) = ^ ~ t«&(«T A + ^1 + tanh^ x( _ R 

(3) 

Note that P 2 (x) + P 2 (—x) = / and therefore P{x) and 
P(—x) describe a measurement. If x = e, where |e| <C 1, 
the measurement is weak. Consider the effect of the oper- 
ators P(x) on a pure state The state can be written 
as \ip) = Pi\ip) + P 2 \ijj) = y/pïtyi) + V^l^), where 
1^1,2) = Pi,2\ip) I '^/Ph2 are the two possible outcomes of 
the projective measurement and pi.2 = (ip\Pi t 2\ip) are 
the corresponding probabilities. If x is positive (nega- 
tive), the operator P{x) increases (decreases) the ratio 
\fp~il \fPx °f the l^) and components of the state. 
By applying the same operator Pis) many times in a row 
for some fixed e, the ratio can be made arbitrarily large 
or small depending on the sign of e, and hence the state 
can be transformed arbitrarily close to |?/>i) or \ip2)- The 
ratio of the p\ and P2 is the only parameter needed to 
describe the state, since p\+p2 = 1 • 

Also note that P(-x)P(x) = (1 - tanh 2 (a;)) 1/2 Í/2 is 
proportional to the identity. If we apply the same mea- 
surement P(±e) twice and two opposite outcomes occur, 
the system returns to its previous state. Thus we see 
that the transformation of the state under many repeti- 
tions of the measurement -P(±e) follows a random walk 
along a curve \ip(x)) in state space. The position on this 
curve can be parameterized by x — ln Wpi/p2· Then 
\ip(x)) can be written as -y/pi (ce) l^i) + \/P2 (íe) IV'2) , where 
p h2 (x) = (l/2)[l±tanh(x)]. 

The measurement given by the operators P(±e) 
changes x by x — > x ± e, with probabilities p±(x) — 
(1 ±t&iíh(e)(pi(x) — p2(x)))/2. We continue this random 
walk until \x\ > X, for some X which is sufhcicntly large 
that \ip(X)) ?» \ipi) and \ip(—X)) w \ip2) to whatever pre- 
cision we desire. What are the respective probabilities of 
these two outcomes? 

Define p(x) to be the probability that the walk will 
end at X (rather than —X) given that it began at x. 
This must satisfy p(x) = p + (x)p{x + e) + p-(x)p(x — e). 
Substituting our expressions for the probabilities, this 
becomes 

p{x) = (p{x + e)+p(x-e))/2 (4) 
+ tanh(e) ta,nh(x)(p(x + e) — p(x — e))/2. 

If we go to the infinitesimal limit e — > dx, this becomes 



a continuous differential equation 
d 2 p , . . dp 

^ + 2tanh(*)J;=0, (5) 

with boundary conditions p(X) = 1, p(— X) = 0. 
The solution to this equation is p(x) — (1/2) [1 + 
tanh(a;)/tanh(X)]. In the limit where X is large, 
tanh(X) — > 1, so p(x) — pi(x). The probabilities of 
the outcomes for the sequence of weak measurements are 
exactly the same as those for a single projective mea- 
surement. Note that this is also true for a walk with a 
step size that is not infinitesimal, since the solution p(x) 
satisfies for an arbitrarily large e. 

Alternatively, instead of looking at the state of the 
system during the process, we could look at an operator 
that effectively describes the system's transformation to 
the current state. This has the advantage that it is state- 
independent, and will lead the way to decompositions of 
generalized measurements; it also becomes obvious that 
the procedure works for mixed states, too. 

We think of the measurement process as a random 
walk along a curve P{x) in operator space, given by 
Eq. ©, which satisfies P(0) = J/-/2, hm P(x) = A, 

x — > — 00 

lim P(x) — P2- It can be verified that P(x)P(y) oc 

x — >oo 

P(x + y), where the constant of proportionality is 
(cosh(a; + y)/2 cosh(x) cosh(y)) 1 / 2 . Due to normalization 
of the state, operators which differ by an overall factor 
are equivalent in their effects on the state. Thus, the ran- 
dom walk driven by weak measurement operators P(±e) 
has a step size |e|. 

Next we consider measurements where the measure- 
ment operators M\ and M2 are positive but not projec- 
tors. We use the well known fact that a generalized mea- 
surement can be implemented as joint unitary operation 
on the system and an ancilla, followed by a projective 
measurement on the ancilla [4|. (One can think of this 
as an indirect measurement; one lets the system interact 
with the ancilla, and then measures the ancilla.) Later 
we will show that the ancilla is not needed. We con- 
sider two-outcome measurements and two-level ancillas. 
In this case M\ and M2 commute, and hence can be si- 
multaneously diagonalized. 

Let the system and ancilla initially be in a state p <g> 
|0)(0|. Consider the unitary operation 

Ü(0) = M x ® Z + M 2 <8> X, (6) 

where X and Z are Pauli matrices acting on the ancilla 
bit. By applying Í7(0) to the extended system we trans- 
form it to: 

Í7(0)(p<8|0}<0|)l7t(0) = M xP Mx® |0> <0| 

+M lP M 2 ® |0)(1| 
+M 2 pMt ® |1)(0| 
+M 2 pM 2 ® (7) 
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Then a projective measurement on the ancilla in the com- 
putational basis would yield one of the possible gener- 
alized measurement outcomes for the system. We can 
perform the projective measurement on the ancilla as a 
sequence of weak measurements by the procedure we de- 
scribed earlier. We will then prové that for this process, 
there exists a corresponding sequence of generalized mea- 
surements with the same effect acting solely on the sys- 
tem. To prové this, we first show that at any stage of 
the measurement process, the state of the extended sys- 
tem can be transformed into the form p(x) ® |0)(0| by a 
unitary operation which does not depend on the state. 

The net effect of the joint unitary operation U(0), fol- 
lowed by the effective measurement operator on the an- 
cilla, can be written in a block form in the computational 
basis of the ancilla: 



M(x) = (ï®P(x))Ü(0) 



Mi 



1 — tanh(:c) 
2 



l+tanh(:c) jÇ^ 



1 — tanh(ic) ji'V 

5 Jyi2 



l+tanh(a;) 
2 



Mi 



(8) 



If the current state M{x){p ® |0)(0|)M^ can be trans- 
formed to p(x)® |0)(0| by a unitary operator U(x) which 

is independent of p, then the lower left block of U (x)M(x) 
should vanish. We look for such a unitary operator in 
block form, with each block being Hermitian and diago- 
nal in the same basis as M and M%. One solution is: 



U(x) 



A(x) È(x) 

è{x) -À(x) 



(9) 



where 



À(x) = y/ï -tanh(x)Mi(Í + tanh(.x)(M 2 2 - M 2 ))~^, 

(10) 

B(x) = v 7 ! + tanh(x)M 2 (í + tanh(a;)(AÍ 2 2 - M 2 ))~i 

(11) 

(Since Mf + M| = J, the operator (I + tanh(a:)(Mf - 
M 2 ))~2 always exists.) Note that U(x) is Hermitian, so 
U (x) = W(x) is its own inverse, and at x = it reduces 
to the operator ((6]). 

After every measurement on the ancilla, depending on 
the value of x, we apply the operation U(x). Then, be- 
fore the next measurement, we apply its inverse U'(x) = 
U(x). By doing this, we can think of the procedure 
as a sequence of generalized measurements on the ex- 
tended system that transform it between states of the 
form p(x) <g> |0)(0| (a generalized measurement preceded 
by a unitary operation and followed by a unitary op- 
eration dependent on the outcome is again a general- 
ized measurement). The measurement operators are now 

M(x,±s) = Ü(x±e)(ï<E)P(±e))Ò(x), and have the form 



M(x,±s) 



M(x,±e) N(x,±e) 
6 Ò{x,±e) 



(12) 



Here M, N, Ò are operators acting on the system. Upon 
measurement, the state of the extended system is trans- 
formed 



p(x) ® 
with probability 



M(x,±s)p(x)MÏ(x,±s) 
p(x,±e) 



O|0)<0|, (13) 



p{x,±e) = Tr{M(x, ±e)p{x)M^ [x,±e)\ . (14) 

By imposing M + (x, e)M(x, e) + M(x, -e)^M(x, -e) = I, 
we obtain that 

M f (x, e)M(x, s) + M 1 (a:, -s)M(x, -e) = Í, (15) 

where the operators in the last equation acts on the sys- 
tem space alone. Therefore, the same transformations 
that the system undergoes during this procedure can be 
achieved by the measurements M(x. ±e) acting solely on 
the system. Depending on the current value of x, we per- 
form the measurement M(x, ±e). Due to the one-to-one 
correspondence with the random walk for the projective 
measurement on the ancilla, this procedure also follows 
a random walk with a step size |e|. It is easy to see that 
if the measurements on the ancilla are weak, the cor- 
responding measurements on the system are also weak. 
Therefore we have shown that every measurement with 
positive operators Mi and M 2 , can be implemented as a 
sequence of weak measurements. This is the main result 
of this paper. From the construction above, one can find 
the explicit form of the weak measurement operators: 



M(x,e) 



1 — tanh(e) 



A(x)À(x + e) 



1 + tanh(e) 



È(x)È(x + e). (16) 



Note that this procedure works even if the step of the 
random walk is not small, since P(x)P{y) cx P(x + y) for 
arbitrary vàlues of x and y. So it is not surprising that 
the effective operator which gives the state of the system 
at the point x is 



M(0,x) = /-^ Miií*) 



1 + tanh(íc) 



M 2 B(x), 



(17) 



where M(x,y) is defincd by (TÏÏÏ)) . 

Finally, consider the most general type of two-outcome 
generalized measurement, with the only restriction being 
Mi Mi + M 2 M2 = By polar decomposition the mea- 
surement operators can be written 



Mi, 2 = ViWMi T i2 Mi, 2 , 



(18) 
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where Vi. 2 are appropriate unitary operators. One can 
think of these unitàries as causing an additional dis- 
turbance to the state of the system, in addition to 
the reduction due to the measurement. The operators 
(M{ 2 Mi i 2) 1 / 2 are positive, and thcy form a measure- 
ment. We could then measure Mi and M 2 by first mea- 
suring these positive operators by a sequence of weak 
measurements, and then performing either V\ or V 2 , de- 
pcnding on the outcome. 

However, we can also decompose this measurement 
directly into a sequence of weak measurements. Let 
the weak measurement operators for {M\ 2 Mi^ 2 ) 1 / 2 be 
M p (x, ±e). Let V{x) be any continuous unitary operator 
function satisfying 1^(0) = I and V(±x) — > Vi t 2 as x — » 

00. We then define M{x,y) = V(x + y)M p (x,y)V" i (x). 
By construction M(x, ±y) are measurement operators. 
Since V(x) is continuous, if y = e, where e <C 1, the 
measurements are weak. The measurement procedure is 
analogous to the previous cases and follows a random 
walk along the curve M(0,x) = V(x)M p (0,x). 

In summary, we have shown that for every two- 
outcome measurement described by operators M% and 
M 2 acting on a Hilbert space of dimension d, there exists 
a continuous two-parameter family of operators M(x,y) 
over the same Hilbert space with the following properties: 
(1) M(x,0) = 7/V2, (2) M(0,x) -> Mi as x -> -00, (3) 
M(0,x) -> M 2 as x -> +00, (4) M(x + y,z)M{x,y) oc 
M(x,z+y), (5) M^x,y)M{x,y)+M^ {x,-y)M{x, -y) 

1 . We have presented an explícit solution for M(x, y) in 
terms of Mi and M 2 . The measurement is implemented 
as a random walk on the curve M(0, x) by consecutive ap- 
plication of the measurements M(x,±e), which depend 
on the current value of the parameter x. In the case 
where |e| -C 1, the measurements driving the random 
walk are weak. Since any measurement can be decom- 
posed into two-outcome measurements, weak measure- 
ments are universal. 

It is obvious from the form of the operators HU) that 
if a measurement is local — the measurement operators 
Mi = Mi ® I act as the identity on all except one sub- 
system of a composite system — it can be implemented as 
sequence of weak local measurements. This result may 
be useful for the study of LOCC (Local Operations with 
Classical Communication) . 

For example, a very useful concept in the theory of en- 
tanglement is the entanglement monotone a function 
of the state that is non-increasing on average under local 
operations. For pure states the operations are unitàries 
and generalized measurements. Since all unitàries can 
be broken into a series of infinitesimal steps and (as we 
have shown) all measurements can be decomposed into 
weak measurements, it sufhces to look at the behavior of 
a prospective monotone under small changes in the state. 
We can thus derive differential conditions for a function 
to be an entanglement monotone. This was one of the 



main motivations for this work. 

Moreover, we can think of the set LOCC itself (or at 
least that subset which preserves pure states) as being 
generated by infinitesimal local operations. This gives 
another way of thinking about entanglement protocols, 
somewhat analogous to studying a Lie àlgebra by exam- 
ining the behavior of its generators. These tòpics will be 
the subject of a follow-up paper [ïo| . 

The connection between weak measurements and 
quantum trajectòries is also an interesting question. 
Quantum trajectòries can be thought of as continuous 
measurements; these generalized measurements should 
therefore correspond to continuous measurements where 
the type of measurement is also continuously adjusted. 
This might be experimentally feasible for some quantum 
optical or atòmic systems, with possible application to 
experiments in quantum control. No doubt the decom- 
position into weak measurements will have many other 
applications; it adds yet another tool to the arsenal of 
the quantum information theorist. 
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